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ABSTRACT 

With the assumption of radial motion and uniform longitudinal distribution 
of coronal mass ejections (CMEs), we propose a method to eliminate projection 
effects from the apparent observed CME latitude distribution. This method has 
been applied to SOHO LASCO data from 1996 January to 2006 December. As 
a result, we find that the real CME latitude distribution had the following char- 
acteristics: (1) High-latitude CMEs {9 > 60° where 9 is the latitude) constituted 
3% of all CMEs and mainly occurred during the time when the polar magnetic 
fields reversed sign. The latitudinal drift of the high-latitude CMEs was corre- 
lated with that of the heliospheric current sheet. (2) 4% of all CMEs occurred in 
the range 45° < 9 < 60°. These mid-latitude CMEs occurred primarily in 2000, 
near the middle of 2002 and in 2005, respectively, forming a prominent three- 
peak structure; (3) The highest occurrence probability of low-latitude {9 < 45°) 
CMEs was at the minimum and during the declining phase of the solar cycle. 
However, the highest occurrence rate of low-latitude CMEs was at the maximum 
and during the declining phase of the solar cycle. The latitudinal evolution of 
low-latitude CMEs did not follow the Sporer sunspot law, which suggests that 
many CMEs originated outside of active regions. 

Subject headings: Sun: coronal mass ejections (CMEs) 



1. Introduction 

Coronal mass ejections (CMEs) are observed with white-light coronagraphs as signifi- 
cant changes in coronal structure moving outward (Skirgiello, 2003). They are believed to 
be the prime driver of most space weather events, e.g., interplanetary shocks, high-energy 
particles and geomagnetic storms (Gosling, 1993). Recently many studies have been carried 



- 2 - 



out on the properties of CME source regions in order to gain insight into the initiation mech- 
anism of such phenomenon. Subramanian and Dere (2001) found that 41% of CME-related 
transients observed are associated with active regions and have no filament eruptions, 44% 
are associated with eruptions of filaments embedded in active regions, and 15% are associated 
with eruptions of filaments outside active regions. Zhou et al. (2003) found that 88% of the 
earth-directed halo CMEs are associated with fiares and 94% are associated with eruptive 
filaments. With regard to the locations of CME source regions, they found that there are 
79% CMEs coming from active regions and 21% originating outside active regions. In this 
letter we concentrate on the latitudinal variations of CME source regions also with a view 
of providing clues about the CME triggering processes. Skirgiello (2003) proposed a mathe- 
matical tool using the apparent latitude distribution to deduce the real latitude distribution. 
By referring to this idea, a much simpler and easier new method with the same function is 
designed in order to obtain the evolution of real CME latitudes during solar cycle 23. 



2. Data 

CMEs are selected from the catalog on the website of http://cdaw.gsfc.nasa.gov/CMEJist 
in the interval 1996 January to 2006 December, as observed by SOHO/LASCO. The SOHO/LASCO 
coronagraph images solar corona continuously since 1996 covering a field of view starting from 
about 1.5 Rs to 32 Rs (Gopalswamy et al. 2003, Rs is the sun's radius). By performing 
preliminary estimation, St. Cyr et al. (2000) found that its detection of CMEs possessed 
an efficiency not less than 95%. The CMEs are observed on the plane of the sky, therefore 
all apparent spatial parameters are a projection of the real values onto that plane. The ap- 
parent latitude 6 comes from the central position angle (CPA), which is defined and listed 
in the CME catalog by Yashiro et al. When 0° < CPA < 180°, 5 = 90° - CPA and when 
180° < CPA < 360°, S = CPA - 270°. The complete halo events with an angular width of 
360° are excluded from the statistics because of difficulty in determining S. 

In consideration of the CME number and the temporal precision, we divide the whole 
interval into 22 parts, as shown in Figure 1. If the time of a part is set very short, the CME 
number will get too small to be able to represent the real distribution oi S; if very long, the 
temporal precision will then get so bad that it is difficult to show the solar cycle variation 
clearly. For the time interval of every part, we compute the occurrence percentages of CMEs 
corresponding to d = 0°, d = 1°, . . ., and S = 90°, respectively. The bin width of S is 1° by 
following that of the CPA. Hence we can define a vector 

S — [so, Si, . . . , Si, . . . , Sgo]'^ (1) 

to give the distribution of 5, where Sj means the percentage at 5 — i°. In the present study. 
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5 is smoothed with a window of 10°, see the thin sohd hnes in Figure 1. Furthermore, the 
distributions for the northern [5 > 0) and southern (S < 0) hemisphere are summed without 
the consideration of the north-south asymmetry of solar activities, e.g., the sunspots (Ohvcr 

6 Ballester, 1994), the filaments (Hansen & Hansen, 1975), the X-ray flares (Li et al. 1998) 
and the high-latitude CMEs (Gopalswamy et al. 2003). 

3. Method 

The proposed method is based on the assumption of radial motion and uniform longi- 
tudinal distribution of CMEs. As shown in Figure 2, the point E indicates the CME source 
region, OE is the direction of CME motion, OH is OE^s projection onto the plane of the 
sky. Thereupon ZHOY, ZLOY, and ZACE are, namely, the apparent latitude d, the real 
latitude 6 and the longitude (/?. By a simple deduction (|Cif| = \CE\sm(p = \0C\/ tan5, 
\CE\ = RgCosO, \0C\ = Rs sinO) we have 

tan 9 — sin (p tan S. (2) 

Set the discrete variable 6 = [0°, 1°, 2°, . . . , 90°]. For each 6 fixed, we can obtain 91 
S values by using Equation 2 if given ip = 0°, (p = l°,---,ip = 90°. We let all 6 equal 
to the expression 5-1-0.5° converted to integer type and then we can define a vector A — 
[oo, Oi, . . . , Oj, . . . , Ogo]^, where Oj means the percentage at S — i° in the obtained 91 S values. 
Some special conditions are dealt with as follows. When ^ = 0°, A = [1, 0, . . . , 0]^; when 
e = 90°, A = [0, . . . , 0, 1]^; when ^ = 0° and ^ ^ 0°, S = 90°. 

Combining 91 vectors A we get a 91 x 91 matrix 

M ^[Ao,Ai,...,Aj,...,Ago], (3) 

where the subscript j means the case 9 — j°. Figure 3 shows the contour of matrix M. If we 
use vector R — [ro, ri, . . . , rgo]^ to indicate the real distribution of ^, we have M x R — S 
and 

R = X S. (4) 

4. Results and Discussion 

The final results of R, also smoothed with a window of 10°, are plotted as 22 thick 
solid hnes in Figure 1. In comparison with S, most high- latitude CMEs are substituted with 
events coming from low latitudes. The reason of this phenomenon is that by projection the 
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CMEs remote from solar limb can be observed at any 5 > 9. For a direct demonstration of 
the solar cycle variation of 6, we combine all vectors R in time order and give a filled contour 
plot in Figure 4a. Figure 4b then shows the variation of the actual (not the percentage) 
CME occurrence rate rj = R x n/t, where n is the CME number, t is the time length (in 
year), both given in Figure 1. 

From Figure 4, we find that most high-latitude CMEs occur during the period from 
1999 to the middle of 2001. Gopalswamy et al. (2003) found a general spreading of CME 
latitudes attaining 60° by 1999. The northern high-latitude CMEs became nonexist beyond 
October 2000. However the southern ones continued to occur until the first quarter of 2002. 
Our time interval of most high-latitude CMEs' occurrence is approximately consistent with 
the above Gopalswamy et al.'s. As shown by white dotted lines, we find a poleward motion 
with a speed of 12.5 deg/year and an equatorward motion with a speed of 25.4 deg/year 
taking place during years 1999-2001 and during 2000 to the middle of 2001, respectively. 
For the poleward motion, years 1999 and 2001 are the polarity reversal time at ^ = 60° and 
the epoch when the polar regions' area is of minimum (Song & Wang, 2006). Therefore we 
suggest this motion come from the collision between the unipolar poleward meridional fiows 
(Song & Wang, 2006) and the polar regions of an opposite polarity. Such collision hastens 
the formation of polar crown filaments whose eruptions have a close correspondence with 
CMEs (Gopalswamy et al. 2003). The much lower speed than that of the meridional flows 
(> 20 deg/year) can be explained by the polar regions blocking the way. The equatorward 
motion may have the same mechanism because during its interval, we can see that the tilt 
angle, defined as the maximum extent of the hehospheric current sheet (HCS) , also declines 
sharply with a similar speed (see Figure 8 in Gopalswamy et al. 2003). The high- latitude 
CMEs constitute 5-10% during the time when the polar magnetic fields reverse sign and 3% 
of all CMEs during the whole solar cycle. 

In the range 9 G [45°, 60°], we only find three prominent peaks located in 2000, near 
the middle of 2002 and in 2005, respectively, which seem to have a period of 2.5 years. By 
comparing these peak positions with the variation of sunspot number (see the pink solid 
lines in Figure 4), we find no any close correlations. Therefore this phenomenon is very 
puzzling indeed. Wang et al. (1989) found many poleward surges of alternating polarities 
in the active belts. Do they initiate the middle-latitude CMEs? We have no definite answer 
at present. Such CMEs constitute 6-12% during their main phases and 4% of all CMEs. 

93% CMEs occur in the range 9 < 45°. From Figure 4, we find that the highest 
occurrence probability (R) of these low-latitude CMEs is at the minimum and during the 
declining phase of solar cycle. However the highest occurrence rate (77) is at the maximum and 
during the declining phase of solar cycle. As shown by green solid lines, at the minimum and 
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during the rising phase, CMEs' average latitude has a hnear increase. Then at the maximum 
and during the dechning phase, CMEs' average latitude is in a nearly steady state, which 
keeps around 15° ~ 25°. It is obvious that it does not follow the Sporer sunspot law. This 
means that there are a large quantity of CMEs originating outside active regions. Zhou 
et al. (2006) classified CME-associated large-scale structures into four different categories: 
extended bipolar regions (EBRs), transequatorial magnetic loops, transequatorial filaments 
and long filaments along the boundaries of EBRs (maybe corresponding to the high-latitude 
CMEs). The latitudes of former three categories are higher or lower than active regions, 
which result in a much wider distribution of CME latitudes. 



5. Conclusions 

Using the proposed method, we have studied the solar cycle variation of real CME 
latitudes during cycle 23 and found three main features. (1) The latitudinal drift of the 
high-latitude CMEs {9 > 60°) was correlated with that of the HCS. (2) The mid-latitude 
CMEs {e e [45°, 60°]) occurred primarily in 2000, near the middle of 2002 and in 2005, 
respectively, forming a prominent three-peak structure. (3) The latitudinal evolution of low- 
latitude CMEs {6 < 45°) did not follow the Sporer sunspot law, which suggests that many 
CMEs originated outside of active regions. Wc believe that such features can provide hints 
for understanding the origin of CMEs. At present, Solar Terrestrial Relations Observatory 
(STEREO), a very nice mission of NASA and launched successfully last October, is observing 
solar corona in 3-D instead of on the plane of the sky. With it the finding of CME source 
regions is no longer a difficult thing and then real CME latitude distribution can be better 
determined. 

This CME catalog is generated and maintained at the CDAW Data Center by NASA 
and The Catholic University of America in cooperation with the Naval Research Laboratory. 
SOHO is a project of international cooperation between ESA and NASA. This work is 
jointly supported by National Natural Science Foundation of China (40621003, 40536029, 
and 40604019), the 973 project under grant 2006CB806304, and the CAS International 
Partnership Program for Creative Research Teams. 
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Fig. 1. — The distributions of the apparent (S, thin sohd hne) and real {R, thick sohd hne) 
CME latitudes during 22 time intervals which are smoothed with a window of 10°. The 
dashed lines around them indicate the unsmoothed variations. In the top right corner of 
each plot, the time range and the CMEs number have been shown. 
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Fig. 2. — The spatial sketch of the apparent latitude {S, ZHOY), the real latitude {9, 
ZLOY) and the longitude {ip, /LACE) of a CME event initiating from the point E. OX is 
the direction of the observer, OE is the direction of CME motion, CA/ /OX, CL/ /OY , and 
EH ± CL. 
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Real latitude, 



Fig. 3. — The contour of the transition matrix M defined in section 3. The distributions 

of the apparent {S) and real (i?) CME latitudes are related by the equation S = M x R. 
Bright white in the upper left and down right corners indicates the value of zero only for 
being pleasing to the eye. The bin widths of 9 and S are both 1°. 
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Fig. 4. — The solar cycle variations of the real CME latitude 9 (a) and the CME occurrence 
rate rj (b). The white dotted lines show the poleward and equatorward motions of the high- 
latitude CMEs. The pink solid line indicates the monthly averaged sunspot number during 
each period of time. The green solid line indicates the average latitude of CMEs in the range 
9 < 45°. Bright white in the upper regions indicates the value of zero. 



